home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_gemm.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  42.7 KB  |  1,400 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_gemm (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int transA = 111;
  14.    int transB = 111;
  15.    int M = 1;
  16.    int N = 2;
  17.    int K = 4;
  18.    float alpha = 1.0f;
  19.    float beta = 0.0f;
  20.    float A[] = { 0.199f, 0.237f, 0.456f, 0.377f };
  21.    int lda = 4;
  22.    float B[] = { 0.842f, -0.734f, 0.323f, -0.957f, -0.303f, -0.873f, -0.871f, -0.819f };
  23.    int ldb = 2;
  24.    float C[] = { 0.498f, -0.925f };
  25.    int ldc = 2;
  26.    float C_expected[] = { -0.222426f, -1.07973f };
  27.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  28.    {
  29.      int i;
  30.      for (i = 0; i < 2; i++) {
  31.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1466)");
  32.      }
  33.    };
  34.   };
  35.  
  36.  
  37.   {
  38.    int order = 102;
  39.    int transA = 111;
  40.    int transB = 111;
  41.    int M = 1;
  42.    int N = 2;
  43.    int K = 4;
  44.    float alpha = 1.0f;
  45.    float beta = 0.0f;
  46.    float A[] = { -0.83f, 0.922f, -0.228f, -0.003f };
  47.    int lda = 1;
  48.    float B[] = { 0.072f, 0.345f, 0.944f, -0.39f, -0.577f, 0.656f, -0.693f, -0.453f };
  49.    int ldb = 4;
  50.    float C[] = { 0.583f, 0.522f };
  51.    int ldc = 1;
  52.    float C_expected[] = { 0.044268f, 1.24311f };
  53.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  54.    {
  55.      int i;
  56.      for (i = 0; i < 2; i++) {
  57.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1467)");
  58.      }
  59.    };
  60.   };
  61.  
  62.  
  63.   {
  64.    int order = 101;
  65.    int transA = 111;
  66.    int transB = 112;
  67.    int M = 1;
  68.    int N = 2;
  69.    int K = 4;
  70.    float alpha = 0.1f;
  71.    float beta = 0.1f;
  72.    float A[] = { -0.838f, 0.622f, -0.494f, 0.304f };
  73.    int lda = 4;
  74.    float B[] = { 0.147f, 0.134f, 0.169f, 0.734f, -0.7f, 0.541f, -0.794f, -0.256f };
  75.    int ldb = 4;
  76.    float C[] = { -0.632f, -0.559f };
  77.    int ldc = 2;
  78.    float C_expected[] = { -0.0532188f, 0.0678514f };
  79.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  80.    {
  81.      int i;
  82.      for (i = 0; i < 2; i++) {
  83.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1468)");
  84.      }
  85.    };
  86.   };
  87.  
  88.  
  89.   {
  90.    int order = 102;
  91.    int transA = 111;
  92.    int transB = 112;
  93.    int M = 1;
  94.    int N = 2;
  95.    int K = 4;
  96.    float alpha = 0.1f;
  97.    float beta = 0.1f;
  98.    float A[] = { -0.937f, 0.635f, 0.596f, -0.51f };
  99.    int lda = 1;
  100.    float B[] = { -0.688f, -0.265f, 0.049f, 0.133f, -0.918f, -0.147f, 0.977f, -0.21f };
  101.    int ldb = 2;
  102.    float C[] = { 0.844f, 0.999f };
  103.    int ldc = 1;
  104.    float C_expected[] = { 0.0474373f, 0.135125f };
  105.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  106.    {
  107.      int i;
  108.      for (i = 0; i < 2; i++) {
  109.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1469)");
  110.      }
  111.    };
  112.   };
  113.  
  114.  
  115.   {
  116.    int order = 101;
  117.    int transA = 112;
  118.    int transB = 111;
  119.    int M = 1;
  120.    int N = 2;
  121.    int K = 4;
  122.    float alpha = -0.3f;
  123.    float beta = 0.1f;
  124.    float A[] = { -0.165f, 0.638f, 0.346f, -0.697f };
  125.    int lda = 1;
  126.    float B[] = { 0.499f, -0.73f, 0.262f, 0.759f, 0.664f, 0.997f, -0.702f, -0.839f };
  127.    int ldb = 2;
  128.    float C[] = { 0.17f, 0.425f };
  129.    int ldc = 2;
  130.    float C_expected[] = { -0.224158f, -0.417831f };
  131.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  132.    {
  133.      int i;
  134.      for (i = 0; i < 2; i++) {
  135.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1470)");
  136.      }
  137.    };
  138.   };
  139.  
  140.  
  141.   {
  142.    int order = 102;
  143.    int transA = 112;
  144.    int transB = 111;
  145.    int M = 1;
  146.    int N = 2;
  147.    int K = 4;
  148.    float alpha = -0.3f;
  149.    float beta = 0.1f;
  150.    float A[] = { -0.603f, -0.714f, -0.893f, 0.046f };
  151.    int lda = 4;
  152.    float B[] = { 0.859f, -0.694f, -0.868f, -0.98f, -0.103f, 0.567f, -0.277f, -0.734f };
  153.    int ldb = 4;
  154.    float C[] = { 0.517f, -0.622f };
  155.    int ldc = 1;
  156.    float C_expected[] = { -0.160575f, -0.0234604f };
  157.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  158.    {
  159.      int i;
  160.      for (i = 0; i < 2; i++) {
  161.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1471)");
  162.      }
  163.    };
  164.   };
  165.  
  166.  
  167.   {
  168.    int order = 101;
  169.    int transA = 112;
  170.    int transB = 112;
  171.    int M = 1;
  172.    int N = 2;
  173.    int K = 4;
  174.    float alpha = 0.1f;
  175.    float beta = 1.0f;
  176.    float A[] = { -0.087f, -0.047f, -0.051f, -0.615f };
  177.    int lda = 1;
  178.    float B[] = { -0.722f, -0.077f, 0.563f, 0.501f, 0.855f, 0.605f, 0.556f, -0.627f };
  179.    int ldb = 4;
  180.    float C[] = { -0.181f, -0.89f };
  181.    int ldc = 2;
  182.    float C_expected[] = { -0.208039f, -0.864557f };
  183.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  184.    {
  185.      int i;
  186.      for (i = 0; i < 2; i++) {
  187.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1472)");
  188.      }
  189.    };
  190.   };
  191.  
  192.  
  193.   {
  194.    int order = 102;
  195.    int transA = 112;
  196.    int transB = 112;
  197.    int M = 1;
  198.    int N = 2;
  199.    int K = 4;
  200.    float alpha = 0.1f;
  201.    float beta = 1.0f;
  202.    float A[] = { -0.753f, -0.074f, -0.247f, -0.19f };
  203.    int lda = 4;
  204.    float B[] = { 0.061f, 0.743f, 0.22f, -0.682f, 0.733f, 0.417f, 0.772f, 0.665f };
  205.    int ldb = 2;
  206.    float C[] = { -0.253f, 0.972f };
  207.    int ldc = 1;
  208.    float C_expected[] = { -0.291994f, 0.898164f };
  209.    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  210.    {
  211.      int i;
  212.      for (i = 0; i < 2; i++) {
  213.        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1473)");
  214.      }
  215.    };
  216.   };
  217.  
  218.  
  219.   {
  220.    int order = 101;
  221.    int transA = 111;
  222.    int transB = 111;
  223.    int M = 1;
  224.    int N = 2;
  225.    int K = 4;
  226.    double alpha = 0;
  227.    double beta = 0;
  228.    double A[] = { 0.017, 0.191, 0.863, -0.97 };
  229.    int lda = 4;
  230.    double B[] = { -0.207, -0.916, -0.278, 0.403, 0.885, 0.409, -0.772, -0.27 };
  231.    int ldb = 2;
  232.    double C[] = { -0.274, -0.858 };
  233.    int ldc = 2;
  234.    double C_expected[] = { 0.0, 0.0 };
  235.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  236.    {
  237.      int i;
  238.      for (i = 0; i < 2; i++) {
  239.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1474)");
  240.      }
  241.    };
  242.   };
  243.  
  244.  
  245.   {
  246.    int order = 102;
  247.    int transA = 111;
  248.    int transB = 111;
  249.    int M = 1;
  250.    int N = 2;
  251.    int K = 4;
  252.    double alpha = 0;
  253.    double beta = 0;
  254.    double A[] = { 0.571, 0.081, 0.109, 0.988 };
  255.    int lda = 1;
  256.    double B[] = { -0.048, -0.753, -0.8, -0.89, -0.535, -0.017, -0.018, -0.544 };
  257.    int ldb = 4;
  258.    double C[] = { -0.876, -0.792 };
  259.    int ldc = 1;
  260.    double C_expected[] = { 0.0, 0.0 };
  261.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  262.    {
  263.      int i;
  264.      for (i = 0; i < 2; i++) {
  265.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1475)");
  266.      }
  267.    };
  268.   };
  269.  
  270.  
  271.   {
  272.    int order = 101;
  273.    int transA = 111;
  274.    int transB = 112;
  275.    int M = 1;
  276.    int N = 2;
  277.    int K = 4;
  278.    double alpha = -0.3;
  279.    double beta = 1;
  280.    double A[] = { 0.939, 0.705, 0.977, 0.4 };
  281.    int lda = 4;
  282.    double B[] = { -0.089, -0.822, 0.937, 0.159, 0.789, -0.413, -0.172, 0.88 };
  283.    int ldb = 4;
  284.    double C[] = { -0.619, 0.063 };
  285.    int ldc = 2;
  286.    double C_expected[] = { -0.7137904, -0.1270986 };
  287.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  288.    {
  289.      int i;
  290.      for (i = 0; i < 2; i++) {
  291.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1476)");
  292.      }
  293.    };
  294.   };
  295.  
  296.  
  297.   {
  298.    int order = 102;
  299.    int transA = 111;
  300.    int transB = 112;
  301.    int M = 1;
  302.    int N = 2;
  303.    int K = 4;
  304.    double alpha = -0.3;
  305.    double beta = 1;
  306.    double A[] = { -0.795, 0.81, 0.388, 0.09 };
  307.    int lda = 1;
  308.    double B[] = { -0.847, 0.031, -0.938, 0.09, -0.286, -0.478, -0.981, 0.881 };
  309.    int ldb = 2;
  310.    double C[] = { -0.242, -0.02 };
  311.    int ldc = 1;
  312.    double C_expected[] = { -0.1562981, -0.0026243 };
  313.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  314.    {
  315.      int i;
  316.      for (i = 0; i < 2; i++) {
  317.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1477)");
  318.      }
  319.    };
  320.   };
  321.  
  322.  
  323.   {
  324.    int order = 101;
  325.    int transA = 112;
  326.    int transB = 111;
  327.    int M = 1;
  328.    int N = 2;
  329.    int K = 4;
  330.    double alpha = -1;
  331.    double beta = 0;
  332.    double A[] = { -0.556, 0.532, 0.746, 0.673 };
  333.    int lda = 1;
  334.    double B[] = { -0.525, 0.967, 0.687, -0.024, 0.527, 0.485, 0.109, -0.46 };
  335.    int ldb = 2;
  336.    double C[] = { -0.495, 0.859 };
  337.    int ldc = 2;
  338.    double C_expected[] = { -1.123883, 0.49819 };
  339.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  340.    {
  341.      int i;
  342.      for (i = 0; i < 2; i++) {
  343.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1478)");
  344.      }
  345.    };
  346.   };
  347.  
  348.  
  349.   {
  350.    int order = 102;
  351.    int transA = 112;
  352.    int transB = 111;
  353.    int M = 1;
  354.    int N = 2;
  355.    int K = 4;
  356.    double alpha = -1;
  357.    double beta = 0;
  358.    double A[] = { -0.358, 0.224, -0.941, 0.513 };
  359.    int lda = 4;
  360.    double B[] = { -0.201, -0.159, -0.586, -0.016, -0.324, 0.411, 0.115, -0.229 };
  361.    int ldb = 4;
  362.    double C[] = { 0.558, 0.596 };
  363.    int ldc = 1;
  364.    double C_expected[] = { -0.57956, 0.017636 };
  365.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  366.    {
  367.      int i;
  368.      for (i = 0; i < 2; i++) {
  369.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1479)");
  370.      }
  371.    };
  372.   };
  373.  
  374.  
  375.   {
  376.    int order = 101;
  377.    int transA = 112;
  378.    int transB = 112;
  379.    int M = 1;
  380.    int N = 2;
  381.    int K = 4;
  382.    double alpha = -0.3;
  383.    double beta = 1;
  384.    double A[] = { -0.586, 0.809, 0.709, -0.524 };
  385.    int lda = 1;
  386.    double B[] = { 0.768, 0.7, 0.619, -0.478, -0.129, -0.778, -0.432, 0.454 };
  387.    int ldb = 4;
  388.    double C[] = { 0.042, 0.252 };
  389.    int ldc = 2;
  390.    double C_expected[] = { -0.1996785, 0.5813976 };
  391.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  392.    {
  393.      int i;
  394.      for (i = 0; i < 2; i++) {
  395.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1480)");
  396.      }
  397.    };
  398.   };
  399.  
  400.  
  401.   {
  402.    int order = 102;
  403.    int transA = 112;
  404.    int transB = 112;
  405.    int M = 1;
  406.    int N = 2;
  407.    int K = 4;
  408.    double alpha = -0.3;
  409.    double beta = 1;
  410.    double A[] = { -0.164, 0.522, 0.948, -0.624 };
  411.    int lda = 4;
  412.    double B[] = { -0.142, 0.778, 0.359, 0.622, -0.637, -0.757, -0.282, -0.805 };
  413.    int ldb = 2;
  414.    double C[] = { -0.09, 0.183 };
  415.    int ldc = 1;
  416.    double C_expected[] = { -0.0248334, 0.1884672 };
  417.    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  418.    {
  419.      int i;
  420.      for (i = 0; i < 2; i++) {
  421.        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1481)");
  422.      }
  423.    };
  424.   };
  425.  
  426.  
  427.   {
  428.    int order = 101;
  429.    int transA = 111;
  430.    int transB = 111;
  431.    int M = 1;
  432.    int N = 2;
  433.    int K = 4;
  434.    float alpha[2] = {0.0f, 1.0f};
  435.    float beta[2] = {0.0f, 0.0f};
  436.    float A[] = { -0.082f, -0.281f, -0.096f, 0.913f, 0.974f, -0.706f, -0.773f, 0.522f };
  437.    int lda = 4;
  438.    float B[] = { 0.745f, -0.664f, 0.352f, -0.733f, 0.304f, -0.555f, -0.493f, -0.089f, 0.188f, 0.631f, 0.235f, 0.152f, -0.299f, -0.731f, -0.686f, -0.332f };
  439.    int ldb = 2;
  440.    float C[] = { -0.179f, -0.284f, -0.996f, -0.414f };
  441.    int ldc = 2;
  442.    float C_expected[] = { -1.06679f, 1.47116f, 0.599689f, 0.933532f };
  443.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  444.    {
  445.      int i;
  446.      for (i = 0; i < 2; i++) {
  447.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1482) real");
  448.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1482) imag");
  449.      };
  450.    };
  451.   };
  452.  
  453.  
  454.   {
  455.    int order = 102;
  456.    int transA = 111;
  457.    int transB = 111;
  458.    int M = 1;
  459.    int N = 2;
  460.    int K = 4;
  461.    float alpha[2] = {0.0f, 1.0f};
  462.    float beta[2] = {0.0f, 0.0f};
  463.    float A[] = { 0.044f, -0.33f, 0.279f, 0.712f, -0.363f, -0.788f, -0.768f, -0.551f };
  464.    int lda = 1;
  465.    float B[] = { 0.138f, 0.927f, -0.178f, -0.864f, 0.888f, 0.844f, -0.199f, 0.706f, -0.034f, 0.483f, 0.499f, 0.664f, 0.648f, 0.324f, 0.97f, 0.609f };
  466.    int ldb = 4;
  467.    float C[] = { -0.129f, 0.842f, 0.214f, -0.626f };
  468.    int ldc = 1;
  469.    float C_expected[] = { 1.81122f, 1.76205f, 1.0574f, -0.564966f };
  470.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  471.    {
  472.      int i;
  473.      for (i = 0; i < 2; i++) {
  474.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1483) real");
  475.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1483) imag");
  476.      };
  477.    };
  478.   };
  479.  
  480.  
  481.   {
  482.    int order = 101;
  483.    int transA = 111;
  484.    int transB = 112;
  485.    int M = 1;
  486.    int N = 2;
  487.    int K = 4;
  488.    float alpha[2] = {0.0f, 0.0f};
  489.    float beta[2] = {-1.0f, 0.0f};
  490.    float A[] = { 0.812f, -0.471f, 0.241f, 0.795f, 0.439f, 0.131f, -0.636f, 0.531f };
  491.    int lda = 4;
  492.    float B[] = { 0.062f, 0.807f, 0.873f, 0.372f, 0.239f, 0.804f, 0.537f, -0.954f, -0.396f, 0.838f, 0.081f, 0.15f, 0.489f, -0.438f, 0.165f, 0.429f };
  493.    int ldb = 4;
  494.    float C[] = { 0.868f, 0.329f, -0.509f, 0.724f };
  495.    int ldc = 2;
  496.    float C_expected[] = { -0.868f, -0.329f, 0.509f, -0.724f };
  497.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  498.    {
  499.      int i;
  500.      for (i = 0; i < 2; i++) {
  501.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1484) real");
  502.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1484) imag");
  503.      };
  504.    };
  505.   };
  506.  
  507.  
  508.   {
  509.    int order = 102;
  510.    int transA = 111;
  511.    int transB = 112;
  512.    int M = 1;
  513.    int N = 2;
  514.    int K = 4;
  515.    float alpha[2] = {0.0f, 0.0f};
  516.    float beta[2] = {-1.0f, 0.0f};
  517.    float A[] = { 0.832f, 0.198f, 0.794f, -0.522f, -0.319f, 0.578f, 0.332f, 0.746f };
  518.    int lda = 1;
  519.    float B[] = { -0.361f, 0.187f, -0.163f, -0.781f, 0.536f, 0.888f, -0.969f, 0.899f, 0.961f, -0.583f, 0.753f, 0.29f, -0.997f, 0.729f, -0.352f, -0.2f };
  520.    int ldb = 2;
  521.    float C[] = { 0.864f, 0.735f, -0.074f, -0.228f };
  522.    int ldc = 1;
  523.    float C_expected[] = { -0.864f, -0.735f, 0.074f, 0.228f };
  524.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  525.    {
  526.      int i;
  527.      for (i = 0; i < 2; i++) {
  528.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1485) real");
  529.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1485) imag");
  530.      };
  531.    };
  532.   };
  533.  
  534.  
  535.   {
  536.    int order = 101;
  537.    int transA = 111;
  538.    int transB = 113;
  539.    int M = 1;
  540.    int N = 2;
  541.    int K = 4;
  542.    float alpha[2] = {0.0f, 1.0f};
  543.    float beta[2] = {0.0f, 0.1f};
  544.    float A[] = { 0.149f, 0.187f, 0.263f, -0.715f, -0.882f, -0.907f, 0.87f, -0.527f };
  545.    int lda = 4;
  546.    float B[] = { -0.915f, -0.249f, -0.986f, -0.799f, -0.136f, 0.712f, 0.964f, 0.799f, -0.569f, 0.686f, 0.603f, 0.758f, 0.161f, -0.698f, -0.263f, -0.256f };
  547.    int ldb = 4;
  548.    float C[] = { 0.622f, -0.824f, -0.482f, -0.161f };
  549.    int ldc = 2;
  550.    float C_expected[] = { -0.246901f, 0.083044f, 1.25556f, 0.009106f };
  551.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  552.    {
  553.      int i;
  554.      for (i = 0; i < 2; i++) {
  555.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1486) real");
  556.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1486) imag");
  557.      };
  558.    };
  559.   };
  560.  
  561.  
  562.   {
  563.    int order = 102;
  564.    int transA = 111;
  565.    int transB = 113;
  566.    int M = 1;
  567.    int N = 2;
  568.    int K = 4;
  569.    float alpha[2] = {0.0f, 1.0f};
  570.    float beta[2] = {0.0f, 0.1f};
  571.    float A[] = { 0.963f, -0.943f, -0.734f, -0.253f, 0.832f, 0.545f, -0.815f, -0.434f };
  572.    int lda = 1;
  573.    float B[] = { 0.23f, -0.211f, 0.906f, 0.232f, -0.339f, 0.597f, -0.919f, 0.793f, 0.535f, 0.526f, 0.119f, 0.053f, 0.751f, 0.044f, 0.752f, -0.469f };
  574.    int ldb = 2;
  575.    float C[] = { 0.483f, -0.266f, -0.224f, -0.692f };
  576.    int ldc = 1;
  577.    float C_expected[] = { -0.047537f, 0.667177f, 1.02025f, 0.823778f };
  578.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  579.    {
  580.      int i;
  581.      for (i = 0; i < 2; i++) {
  582.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1487) real");
  583.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1487) imag");
  584.      };
  585.    };
  586.   };
  587.  
  588.  
  589.   {
  590.    int order = 101;
  591.    int transA = 112;
  592.    int transB = 111;
  593.    int M = 1;
  594.    int N = 2;
  595.    int K = 4;
  596.    float alpha[2] = {-0.3f, 0.1f};
  597.    float beta[2] = {-1.0f, 0.0f};
  598.    float A[] = { -0.657f, -0.497f, -0.293f, -0.168f, -0.943f, -0.181f, 0.569f, 0.91f };
  599.    int lda = 1;
  600.    float B[] = { -0.047f, 0.796f, -0.913f, 0.998f, 0.365f, 0.467f, -0.627f, -0.523f, 0.885f, 0.234f, -0.494f, 0.071f, -0.361f, -0.154f, -0.055f, -0.32f };
  601.    int ldb = 2;
  602.    float C[] = { 0.956f, 0.268f, 0.152f, 0.717f };
  603.    int ldc = 2;
  604.    float C_expected[] = { -0.668685f, 0.134477f, -0.715786f, -0.478065f };
  605.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  606.    {
  607.      int i;
  608.      for (i = 0; i < 2; i++) {
  609.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1488) real");
  610.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1488) imag");
  611.      };
  612.    };
  613.   };
  614.  
  615.  
  616.   {
  617.    int order = 102;
  618.    int transA = 112;
  619.    int transB = 111;
  620.    int M = 1;
  621.    int N = 2;
  622.    int K = 4;
  623.    float alpha[2] = {-0.3f, 0.1f};
  624.    float beta[2] = {-1.0f, 0.0f};
  625.    float A[] = { 0.394f, -0.482f, 0.631f, -0.833f, 0.221f, 0.672f, 0.2f, 0.967f };
  626.    int lda = 4;
  627.    float B[] = { 0.708f, 0.695f, 0.111f, -0.912f, 0.376f, 0.606f, -0.997f, -0.741f, 0.349f, 0.543f, 0.372f, -0.563f, 0.129f, -0.295f, -0.672f, -0.95f };
  628.    int ldb = 4;
  629.    float C[] = { 0.436f, 0.752f, 0.074f, 0.209f };
  630.    int ldc = 1;
  631.    float C_expected[] = { -0.325083f, -0.301952f, -0.283022f, 0.339919f };
  632.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  633.    {
  634.      int i;
  635.      for (i = 0; i < 2; i++) {
  636.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1489) real");
  637.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1489) imag");
  638.      };
  639.    };
  640.   };
  641.  
  642.  
  643.   {
  644.    int order = 101;
  645.    int transA = 112;
  646.    int transB = 112;
  647.    int M = 1;
  648.    int N = 2;
  649.    int K = 4;
  650.    float alpha[2] = {1.0f, 0.0f};
  651.    float beta[2] = {-0.3f, 0.1f};
  652.    float A[] = { 0.827f, -0.862f, 0.373f, -0.265f, -0.9f, 0.892f, -0.319f, 0.151f };
  653.    int lda = 1;
  654.    float B[] = { 0.603f, 0.816f, -0.511f, 0.831f, -0.36f, -0.954f, -0.978f, 0.485f, 0.675f, 0.186f, 0.463f, 0.144f, 0.851f, -0.458f, 0.766f, -0.213f };
  655.    int ldb = 4;
  656.    float C[] = { -0.335f, 0.333f, -0.4f, 0.422f };
  657.    int ldc = 2;
  658.    float C_expected[] = { 2.7126f, 0.702111f, 0.437661f, 0.691294f };
  659.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  660.    {
  661.      int i;
  662.      for (i = 0; i < 2; i++) {
  663.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1490) real");
  664.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1490) imag");
  665.      };
  666.    };
  667.   };
  668.  
  669.  
  670.   {
  671.    int order = 102;
  672.    int transA = 112;
  673.    int transB = 112;
  674.    int M = 1;
  675.    int N = 2;
  676.    int K = 4;
  677.    float alpha[2] = {1.0f, 0.0f};
  678.    float beta[2] = {-0.3f, 0.1f};
  679.    float A[] = { 0.966f, 0.476f, -0.013f, -0.655f, 0.773f, -0.543f, -0.231f, -0.353f };
  680.    int lda = 4;
  681.    float B[] = { -0.684f, 0.144f, 0.018f, -0.77f, -0.688f, 0.909f, -0.094f, -0.938f, -0.757f, 0.574f, -0.479f, 0.473f, 0.0f, 0.064f, -0.168f, 0.858f };
  682.    int ldb = 2;
  683.    float C[] = { -0.912f, 0.54f, 0.756f, 0.024f };
  684.    int ldc = 1;
  685.    float C_expected[] = { -0.156236f, 0.839112f, -0.230206f, -0.106256f };
  686.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  687.    {
  688.      int i;
  689.      for (i = 0; i < 2; i++) {
  690.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1491) real");
  691.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1491) imag");
  692.      };
  693.    };
  694.   };
  695.  
  696.  
  697.   {
  698.    int order = 101;
  699.    int transA = 112;
  700.    int transB = 113;
  701.    int M = 1;
  702.    int N = 2;
  703.    int K = 4;
  704.    float alpha[2] = {0.0f, 0.0f};
  705.    float beta[2] = {1.0f, 0.0f};
  706.    float A[] = { 0.66f, -0.113f, -0.663f, -0.856f, 0.614f, -0.344f, -0.964f, -0.532f };
  707.    int lda = 1;
  708.    float B[] = { -0.606f, -0.965f, -0.279f, -0.312f, 0.63f, 0.967f, 0.041f, -0.557f, 0.663f, 0.619f, -0.134f, 0.261f, -0.388f, 0.525f, 0.222f, 0.538f };
  709.    int ldb = 4;
  710.    float C[] = { 0.114f, -0.376f, -0.851f, -0.682f };
  711.    int ldc = 2;
  712.    float C_expected[] = { 0.114f, -0.376f, -0.851f, -0.682f };
  713.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  714.    {
  715.      int i;
  716.      for (i = 0; i < 2; i++) {
  717.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1492) real");
  718.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1492) imag");
  719.      };
  720.    };
  721.   };
  722.  
  723.  
  724.   {
  725.    int order = 102;
  726.    int transA = 112;
  727.    int transB = 113;
  728.    int M = 1;
  729.    int N = 2;
  730.    int K = 4;
  731.    float alpha[2] = {0.0f, 0.0f};
  732.    float beta[2] = {1.0f, 0.0f};
  733.    float A[] = { 0.212f, -0.752f, 0.679f, 0.49f, -0.029f, -0.488f, 0.567f, 0.374f };
  734.    int lda = 4;
  735.    float B[] = { -0.914f, 0.734f, -0.845f, 0.059f, -0.297f, 0.152f, -0.417f, -0.669f, 0.831f, -0.544f, 0.022f, 0.102f, -0.379f, -0.357f, -0.394f, -0.588f };
  736.    int ldb = 2;
  737.    float C[] = { -0.584f, 0.373f, 0.235f, 0.521f };
  738.    int ldc = 1;
  739.    float C_expected[] = { -0.584f, 0.373f, 0.235f, 0.521f };
  740.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  741.    {
  742.      int i;
  743.      for (i = 0; i < 2; i++) {
  744.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1493) real");
  745.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1493) imag");
  746.      };
  747.    };
  748.   };
  749.  
  750.  
  751.   {
  752.    int order = 101;
  753.    int transA = 113;
  754.    int transB = 111;
  755.    int M = 1;
  756.    int N = 2;
  757.    int K = 4;
  758.    float alpha[2] = {0.0f, 0.1f};
  759.    float beta[2] = {-1.0f, 0.0f};
  760.    float A[] = { 0.135f, 0.128f, 0.909f, -0.963f, 0.299f, -0.944f, 0.944f, 0.942f };
  761.    int lda = 1;
  762.    float B[] = { 0.924f, -0.317f, -0.992f, -0.854f, -0.435f, 0.102f, 0.126f, 0.862f, 0.952f, 0.68f, 0.545f, 0.168f, 0.752f, 0.549f, 0.687f, -0.76f };
  763.    int ldb = 2;
  764.    float C[] = { -0.369f, -0.33f, 0.849f, -0.632f };
  765.    int ldc = 2;
  766.    float C_expected[] = { 0.326537f, 0.37603f, -0.86067f, 0.529817f };
  767.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  768.    {
  769.      int i;
  770.      for (i = 0; i < 2; i++) {
  771.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1494) real");
  772.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1494) imag");
  773.      };
  774.    };
  775.   };
  776.  
  777.  
  778.   {
  779.    int order = 102;
  780.    int transA = 113;
  781.    int transB = 111;
  782.    int M = 1;
  783.    int N = 2;
  784.    int K = 4;
  785.    float alpha[2] = {0.0f, 0.1f};
  786.    float beta[2] = {-1.0f, 0.0f};
  787.    float A[] = { 0.061f, -0.271f, -0.043f, -0.023f, 0.694f, 0.333f, 0.733f, -0.967f };
  788.    int lda = 4;
  789.    float B[] = { 0.088f, -0.607f, 0.589f, 0.375f, -0.897f, -0.954f, -0.216f, -0.195f, -0.865f, -0.511f, -0.219f, 0.535f, 0.976f, 0.582f, 0.464f, -0.041f };
  790.    int ldb = 4;
  791.    float C[] = { 0.533f, -0.63f, 0.405f, 0.667f };
  792.    int ldc = 1;
  793.    float C_expected[] = { -0.459906f, 0.552595f, -0.425391f, -0.533626f };
  794.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  795.    {
  796.      int i;
  797.      for (i = 0; i < 2; i++) {
  798.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1495) real");
  799.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1495) imag");
  800.      };
  801.    };
  802.   };
  803.  
  804.  
  805.   {
  806.    int order = 101;
  807.    int transA = 113;
  808.    int transB = 112;
  809.    int M = 1;
  810.    int N = 2;
  811.    int K = 4;
  812.    float alpha[2] = {0.0f, 0.0f};
  813.    float beta[2] = {1.0f, 0.0f};
  814.    float A[] = { -0.676f, -0.116f, 0.707f, -0.256f, -0.893f, -0.966f, 0.159f, -0.246f };
  815.    int lda = 1;
  816.    float B[] = { 0.059f, 0.281f, -0.93f, -0.263f, 0.583f, -0.11f, 0.639f, -0.96f, -0.878f, 0.984f, 0.058f, 0.977f, -0.567f, 0.561f, -0.048f, -0.798f };
  817.    int ldb = 4;
  818.    float C[] = { 0.362f, -0.808f, 0.428f, -0.112f };
  819.    int ldc = 2;
  820.    float C_expected[] = { 0.362f, -0.808f, 0.428f, -0.112f };
  821.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  822.    {
  823.      int i;
  824.      for (i = 0; i < 2; i++) {
  825.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1496) real");
  826.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1496) imag");
  827.      };
  828.    };
  829.   };
  830.  
  831.  
  832.   {
  833.    int order = 102;
  834.    int transA = 113;
  835.    int transB = 112;
  836.    int M = 1;
  837.    int N = 2;
  838.    int K = 4;
  839.    float alpha[2] = {0.0f, 0.0f};
  840.    float beta[2] = {1.0f, 0.0f};
  841.    float A[] = { -0.915f, 0.439f, 0.171f, -0.019f, 0.843f, 0.944f, -0.581f, 0.856f };
  842.    int lda = 4;
  843.    float B[] = { -0.284f, 0.207f, -0.27f, 0.832f, 0.894f, -0.626f, -0.305f, -0.006f, 0.562f, -0.744f, -0.533f, 0.126f, -0.375f, -0.333f, 0.275f, 0.748f };
  844.    int ldb = 2;
  845.    float C[] = { -0.763f, -0.829f, 0.708f, -0.613f };
  846.    int ldc = 1;
  847.    float C_expected[] = { -0.763f, -0.829f, 0.708f, -0.613f };
  848.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  849.    {
  850.      int i;
  851.      for (i = 0; i < 2; i++) {
  852.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1497) real");
  853.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1497) imag");
  854.      };
  855.    };
  856.   };
  857.  
  858.  
  859.   {
  860.    int order = 101;
  861.    int transA = 113;
  862.    int transB = 113;
  863.    int M = 1;
  864.    int N = 2;
  865.    int K = 4;
  866.    float alpha[2] = {0.0f, 0.1f};
  867.    float beta[2] = {0.0f, 1.0f};
  868.    float A[] = { 0.496f, -0.9f, 0.825f, -0.678f, 0.41f, -0.585f, -0.264f, 0.308f };
  869.    int lda = 1;
  870.    float B[] = { 0.907f, 0.972f, -0.724f, 0.745f, -0.601f, 0.589f, 0.759f, -0.521f, -0.161f, -0.321f, 0.341f, -0.981f, -0.378f, -0.671f, -0.314f, -0.878f };
  871.    int ldb = 4;
  872.    float C[] = { -0.293f, 0.07f, 0.087f, -0.542f };
  873.    int ldc = 2;
  874.    float C_expected[] = { 0.10357f, -0.163927f, 0.444626f, -0.0076744f };
  875.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  876.    {
  877.      int i;
  878.      for (i = 0; i < 2; i++) {
  879.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1498) real");
  880.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1498) imag");
  881.      };
  882.    };
  883.   };
  884.  
  885.  
  886.   {
  887.    int order = 102;
  888.    int transA = 113;
  889.    int transB = 113;
  890.    int M = 1;
  891.    int N = 2;
  892.    int K = 4;
  893.    float alpha[2] = {0.0f, 0.1f};
  894.    float beta[2] = {0.0f, 1.0f};
  895.    float A[] = { -0.225f, -0.629f, -0.939f, -0.836f, -0.841f, -0.794f, 0.836f, -0.65f };
  896.    int lda = 4;
  897.    float B[] = { 0.869f, -0.453f, 0.8f, -0.947f, 0.545f, 0.716f, -0.507f, -0.228f, 0.722f, 0.372f, 0.77f, 0.317f, -0.153f, -0.524f, -0.465f, -0.684f };
  898.    int ldb = 2;
  899.    float C[] = { -0.896f, 0.91f, -0.973f, -0.269f };
  900.    int ldc = 1;
  901.    float C_expected[] = { -1.18974f, -1.0134f, 0.189027f, -1.14494f };
  902.    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  903.    {
  904.      int i;
  905.      for (i = 0; i < 2; i++) {
  906.        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1499) real");
  907.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1499) imag");
  908.      };
  909.    };
  910.   };
  911.  
  912.  
  913.   {
  914.    int order = 101;
  915.    int transA = 111;
  916.    int transB = 111;
  917.    int M = 1;
  918.    int N = 2;
  919.    int K = 4;
  920.    double alpha[2] = {1, 0};
  921.    double beta[2] = {-1, 0};
  922.    double A[] = { -0.33, 0.457, 0.428, -0.19, 0.86, -0.53, 0.058, -0.942 };
  923.    int lda = 4;
  924.    double B[] = { 0.434, 0.653, -0.124, 0.191, -0.112, -0.84, -0.72, 0.075, -0.503, -0.109, 0.3, -0.898, 0.489, 0.384, 0.993, -0.804 };
  925.    int ldb = 2;
  926.    double C[] = { -0.792, -0.155, -0.608, -0.243 };
  927.    int ldc = 2;
  928.    double C_expected[] = { 0.042563, -0.465908, -0.649991, -1.621116 };
  929.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  930.    {
  931.      int i;
  932.      for (i = 0; i < 2; i++) {
  933.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1500) real");
  934.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1500) imag");
  935.      };
  936.    };
  937.   };
  938.  
  939.  
  940.   {
  941.    int order = 102;
  942.    int transA = 111;
  943.    int transB = 111;
  944.    int M = 1;
  945.    int N = 2;
  946.    int K = 4;
  947.    double alpha[2] = {1, 0};
  948.    double beta[2] = {-1, 0};
  949.    double A[] = { 0.726, -0.438, -0.23, -0.054, -0.019, 0.902, -0.883, -0.235 };
  950.    int lda = 1;
  951.    double B[] = { 0.159, -0.18, 0.386, -0.167, 0.971, -0.072, 0.87, -0.839, 0.474, 0.956, -0.235, 0.332, 0.826, -0.056, -0.941, 0.01 };
  952.    int ldb = 4;
  953.    double C[] = { -0.799, 0.973, -0.549, -0.177 };
  954.    int ldc = 1;
  955.    double C_expected[] = { -0.181084, 0.257841, 2.251901, 1.558195 };
  956.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  957.    {
  958.      int i;
  959.      for (i = 0; i < 2; i++) {
  960.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1501) real");
  961.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1501) imag");
  962.      };
  963.    };
  964.   };
  965.  
  966.  
  967.   {
  968.    int order = 101;
  969.    int transA = 111;
  970.    int transB = 112;
  971.    int M = 1;
  972.    int N = 2;
  973.    int K = 4;
  974.    double alpha[2] = {0, 0.1};
  975.    double beta[2] = {1, 0};
  976.    double A[] = { 0.109, 0.892, -0.723, 0.793, 0.109, -0.419, -0.534, 0.448 };
  977.    int lda = 4;
  978.    double B[] = { -0.875, -0.31, -0.027, 0.067, 0.274, -0.126, -0.548, 0.497, 0.681, 0.388, 0.909, 0.889, 0.982, -0.074, -0.788, 0.233 };
  979.    int ldb = 4;
  980.    double C[] = { 0.503, 0.067, 0.239, 0.876 };
  981.    int ldc = 2;
  982.    double C_expected[] = { 0.6553584, 0.0864583, 0.2559136, 0.7518389 };
  983.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  984.    {
  985.      int i;
  986.      for (i = 0; i < 2; i++) {
  987.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1502) real");
  988.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1502) imag");
  989.      };
  990.    };
  991.   };
  992.  
  993.  
  994.   {
  995.    int order = 102;
  996.    int transA = 111;
  997.    int transB = 112;
  998.    int M = 1;
  999.    int N = 2;
  1000.    int K = 4;
  1001.    double alpha[2] = {0, 0.1};
  1002.    double beta[2] = {1, 0};
  1003.    double A[] = { 0.334, 0.192, -0.992, -0.168, 0.154, -0.75, -0.797, -0.76 };
  1004.    int lda = 1;
  1005.    double B[] = { -0.82, 0.147, -0.237, 0.68, 0.317, 0.257, -0.406, -0.802, 0.058, 0.012, -0.832, 0.949, -0.263, -0.085, -0.064, 0.492 };
  1006.    int ldb = 2;
  1007.    double C[] = { 0.079, -0.602, -0.392, 0.316 };
  1008.    int ldc = 1;
  1009.    double C_expected[] = { 0.0980569, -0.6430449, -0.539207, 0.4226848 };
  1010.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1011.    {
  1012.      int i;
  1013.      for (i = 0; i < 2; i++) {
  1014.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1503) real");
  1015.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1503) imag");
  1016.      };
  1017.    };
  1018.   };
  1019.  
  1020.  
  1021.   {
  1022.    int order = 101;
  1023.    int transA = 111;
  1024.    int transB = 113;
  1025.    int M = 1;
  1026.    int N = 2;
  1027.    int K = 4;
  1028.    double alpha[2] = {0, 0};
  1029.    double beta[2] = {-1, 0};
  1030.    double A[] = { -0.305, -0.698, -0.072, -0.383, 0.364, -0.656, 0.819, 0.194 };
  1031.    int lda = 4;
  1032.    double B[] = { 0.682, 0.498, -0.389, 0.923, -0.853, -0.558, -0.722, -0.085, -0.27, 0.026, -0.107, -0.036, 0.644, -0.327, -0.894, 0.34 };
  1033.    int ldb = 4;
  1034.    double C[] = { 0.981, -0.336, -0.377, -0.41 };
  1035.    int ldc = 2;
  1036.    double C_expected[] = { -0.981, 0.336, 0.377, 0.41 };
  1037.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1038.    {
  1039.      int i;
  1040.      for (i = 0; i < 2; i++) {
  1041.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1504) real");
  1042.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1504) imag");
  1043.      };
  1044.    };
  1045.   };
  1046.  
  1047.  
  1048.   {
  1049.    int order = 102;
  1050.    int transA = 111;
  1051.    int transB = 113;
  1052.    int M = 1;
  1053.    int N = 2;
  1054.    int K = 4;
  1055.    double alpha[2] = {0, 0};
  1056.    double beta[2] = {-1, 0};
  1057.    double A[] = { -0.306, -0.709, -0.196, 0.285, 0.873, -0.802, 0.715, -0.179 };
  1058.    int lda = 1;
  1059.    double B[] = { 0.028, 0.109, 0.87, -0.446, 0.735, 0.731, 0.021, -0.186, 0.541, 0.97, -0.333, 0.002, -0.089, -0.01, 0.331, 0.851 };
  1060.    int ldb = 2;
  1061.    double C[] = { 0.902, -0.584, -0.695, -0.607 };
  1062.    int ldc = 1;
  1063.    double C_expected[] = { -0.902, 0.584, 0.695, 0.607 };
  1064.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1065.    {
  1066.      int i;
  1067.      for (i = 0; i < 2; i++) {
  1068.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1505) real");
  1069.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1505) imag");
  1070.      };
  1071.    };
  1072.   };
  1073.  
  1074.  
  1075.   {
  1076.    int order = 101;
  1077.    int transA = 112;
  1078.    int transB = 111;
  1079.    int M = 1;
  1080.    int N = 2;
  1081.    int K = 4;
  1082.    double alpha[2] = {-1, 0};
  1083.    double beta[2] = {1, 0};
  1084.    double A[] = { 0.517, -0.136, 0.72, -0.237, 0.121, -0.66, 0.005, 0.759 };
  1085.    int lda = 1;
  1086.    double B[] = { -0.606, 0.049, 0.807, -0.236, -0.258, -0.412, 0.75, -0.659, 0.993, -0.029, -0.968, 0.707, -0.362, -0.005, 0.096, -0.241 };
  1087.    int ldb = 2;
  1088.    double C[] = { 0.63, 0.922, 0.025, -0.535 };
  1089.    int ldc = 2;
  1090.    double C_expected[] = { 1.117044, 1.983417, -1.276831, -0.447092 };
  1091.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1092.    {
  1093.      int i;
  1094.      for (i = 0; i < 2; i++) {
  1095.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1506) real");
  1096.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1506) imag");
  1097.      };
  1098.    };
  1099.   };
  1100.  
  1101.  
  1102.   {
  1103.    int order = 102;
  1104.    int transA = 112;
  1105.    int transB = 111;
  1106.    int M = 1;
  1107.    int N = 2;
  1108.    int K = 4;
  1109.    double alpha[2] = {-1, 0};
  1110.    double beta[2] = {1, 0};
  1111.    double A[] = { 0.064, 0.371, -0.01, -0.262, 0.143, -0.081, 0.1, -0.062 };
  1112.    int lda = 4;
  1113.    double B[] = { -0.749, 0.289, -0.239, -0.226, 0.284, 0.668, 0.305, 0.075, -0.36, 0.166, -0.416, 0.234, -0.267, 0.525, 0.116, -0.561 };
  1114.    int ldb = 4;
  1115.    double C[] = { 0.671, 0.763, 0.444, -0.246 };
  1116.    int ldc = 1;
  1117.    double C_expected[] = { 0.753107, 0.896395, 0.481996, -0.263126 };
  1118.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1119.    {
  1120.      int i;
  1121.      for (i = 0; i < 2; i++) {
  1122.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1507) real");
  1123.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1507) imag");
  1124.      };
  1125.    };
  1126.   };
  1127.  
  1128.  
  1129.   {
  1130.    int order = 101;
  1131.    int transA = 112;
  1132.    int transB = 112;
  1133.    int M = 1;
  1134.    int N = 2;
  1135.    int K = 4;
  1136.    double alpha[2] = {1, 0};
  1137.    double beta[2] = {-0.3, 0.1};
  1138.    double A[] = { -0.956, -0.751, 0.671, -0.633, 0.648, -0.042, 0.948, 0.826 };
  1139.    int lda = 1;
  1140.    double B[] = { 0.921, 0.506, -0.609, 0.817, -0.686, 0.991, 0.616, -0.482, -0.02, -0.34, 0.559, 0.976, 0.431, 0.385, -0.164, -0.778 };
  1141.    int ldb = 4;
  1142.    double C[] = { 0.074, -0.01, 0.165, 0.166 };
  1143.    int ldc = 2;
  1144.    double C_expected[] = { 0.166046, 0.491557, 1.473191, -0.033821 };
  1145.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1146.    {
  1147.      int i;
  1148.      for (i = 0; i < 2; i++) {
  1149.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1508) real");
  1150.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1508) imag");
  1151.      };
  1152.    };
  1153.   };
  1154.  
  1155.  
  1156.   {
  1157.    int order = 102;
  1158.    int transA = 112;
  1159.    int transB = 112;
  1160.    int M = 1;
  1161.    int N = 2;
  1162.    int K = 4;
  1163.    double alpha[2] = {1, 0};
  1164.    double beta[2] = {-0.3, 0.1};
  1165.    double A[] = { -0.698, -0.062, 0.023, 0.704, 0.443, -0.46, 0.541, 0.296 };
  1166.    int lda = 4;
  1167.    double B[] = { 0.787, -0.199, 0.835, -0.276, -0.515, 0.467, -0.76, -0.483, 0.015, -0.394, -0.748, 0.02, 0.573, 0.3, -0.088, -0.238 };
  1168.    int ldb = 2;
  1169.    double C[] = { 0.935, -0.655, -0.797, 0.071 };
  1170.    int ldc = 1;
  1171.    double C_expected[] = { -1.070679, 0.178755, -0.344714, -0.308137 };
  1172.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1173.    {
  1174.      int i;
  1175.      for (i = 0; i < 2; i++) {
  1176.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1509) real");
  1177.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1509) imag");
  1178.      };
  1179.    };
  1180.   };
  1181.  
  1182.  
  1183.   {
  1184.    int order = 101;
  1185.    int transA = 112;
  1186.    int transB = 113;
  1187.    int M = 1;
  1188.    int N = 2;
  1189.    int K = 4;
  1190.    double alpha[2] = {0, 0.1};
  1191.    double beta[2] = {-0.3, 0.1};
  1192.    double A[] = { -0.202, -0.219, 0.741, 0.527, 0.054, 0.16, -0.359, 0.338 };
  1193.    int lda = 1;
  1194.    double B[] = { -0.872, 0.995, 0.722, 0.618, -0.27, 0.939, -0.743, 0.547, -0.864, 0.376, -0.997, -0.63, 0.887, -0.454, 0.436, -0.039 };
  1195.    int ldb = 4;
  1196.    double C[] = { -0.684, 0.463, -0.386, -0.524 };
  1197.    int ldc = 2;
  1198.    double C_expected[] = { 0.1423153, -0.066679, 0.1175618, 0.0012949 };
  1199.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1200.    {
  1201.      int i;
  1202.      for (i = 0; i < 2; i++) {
  1203.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1510) real");
  1204.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1510) imag");
  1205.      };
  1206.    };
  1207.   };
  1208.  
  1209.  
  1210.   {
  1211.    int order = 102;
  1212.    int transA = 112;
  1213.    int transB = 113;
  1214.    int M = 1;
  1215.    int N = 2;
  1216.    int K = 4;
  1217.    double alpha[2] = {0, 0.1};
  1218.    double beta[2] = {-0.3, 0.1};
  1219.    double A[] = { -0.855, -0.173, -0.679, 0.824, 0.469, 0.786, 0.757, -0.109 };
  1220.    int lda = 4;
  1221.    double B[] = { 0.483, -0.888, -0.757, 0.551, -0.81, 0.23, -0.078, 0.725, -0.592, 0.394, 0.884, 0.802, -0.813, -0.016, -0.853, 0.783 };
  1222.    int ldb = 2;
  1223.    double C[] = { 0.181, -0.368, -0.864, -0.784 };
  1224.    int ldc = 1;
  1225.    double C_expected[] = { 0.1728438, 0.1183508, 0.2526999, 0.3004174 };
  1226.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1227.    {
  1228.      int i;
  1229.      for (i = 0; i < 2; i++) {
  1230.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1511) real");
  1231.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1511) imag");
  1232.      };
  1233.    };
  1234.   };
  1235.  
  1236.  
  1237.   {
  1238.    int order = 101;
  1239.    int transA = 113;
  1240.    int transB = 111;
  1241.    int M = 1;
  1242.    int N = 2;
  1243.    int K = 4;
  1244.    double alpha[2] = {-1, 0};
  1245.    double beta[2] = {-0.3, 0.1};
  1246.    double A[] = { 0.446, -0.65, -0.724, 0.014, 0.792, -0.695, -0.81, -0.358 };
  1247.    int lda = 1;
  1248.    double B[] = { -0.08, 0.216, 0.689, 0.699, 0.073, -0.346, 0.821, -0.668, -0.798, 0.869, 0.451, -0.061, -0.41, 0.316, 0.104, -0.514 };
  1249.    int ldb = 2;
  1250.    double C[] = { -0.476, 0.211, -0.912, -0.243 };
  1251.    int ldc = 2;
  1252.    double C_expected[] = { 1.372475, -0.135616, 0.549353, -1.968747 };
  1253.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1254.    {
  1255.      int i;
  1256.      for (i = 0; i < 2; i++) {
  1257.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1512) real");
  1258.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1512) imag");
  1259.      };
  1260.    };
  1261.   };
  1262.  
  1263.  
  1264.   {
  1265.    int order = 102;
  1266.    int transA = 113;
  1267.    int transB = 111;
  1268.    int M = 1;
  1269.    int N = 2;
  1270.    int K = 4;
  1271.    double alpha[2] = {-1, 0};
  1272.    double beta[2] = {-0.3, 0.1};
  1273.    double A[] = { 0.669, 0.046, -0.094, 0.666, 0.23, 0.448, -0.795, -0.142 };
  1274.    int lda = 4;
  1275.    double B[] = { 0.037, -0.154, -0.739, 0.905, 0.793, -0.53, -0.34, 0.428, 0.072, -0.263, -0.603, -0.905, 0.681, -0.083, -0.511, -0.337 };
  1276.    int ldb = 4;
  1277.    double C[] = { 0.247, 0.575, -0.836, -0.883 };
  1278.    int ldc = 1;
  1279.    double C_expected[] = { -0.975939, 0.415528, 0.275533, 0.002716 };
  1280.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1281.    {
  1282.      int i;
  1283.      for (i = 0; i < 2; i++) {
  1284.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1513) real");
  1285.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1513) imag");
  1286.      };
  1287.    };
  1288.   };
  1289.  
  1290.  
  1291.   {
  1292.    int order = 101;
  1293.    int transA = 113;
  1294.    int transB = 112;
  1295.    int M = 1;
  1296.    int N = 2;
  1297.    int K = 4;
  1298.    double alpha[2] = {0, 0};
  1299.    double beta[2] = {-1, 0};
  1300.    double A[] = { 0.369, 0.506, 0.217, -0.739, -0.395, 0.16, -0.329, -0.954 };
  1301.    int lda = 1;
  1302.    double B[] = { -0.622, -0.945, 0.416, -0.884, 0.797, -0.74, 0.519, -0.789, -0.348, 0.563, -0.398, -0.956, 0.227, 0.84, -0.079, 0.847 };
  1303.    int ldb = 4;
  1304.    double C[] = { 0.833, 0.761, 0.074, -0.448 };
  1305.    int ldc = 2;
  1306.    double C_expected[] = { -0.833, -0.761, -0.074, 0.448 };
  1307.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1308.    {
  1309.      int i;
  1310.      for (i = 0; i < 2; i++) {
  1311.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1514) real");
  1312.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1514) imag");
  1313.      };
  1314.    };
  1315.   };
  1316.  
  1317.  
  1318.   {
  1319.    int order = 102;
  1320.    int transA = 113;
  1321.    int transB = 112;
  1322.    int M = 1;
  1323.    int N = 2;
  1324.    int K = 4;
  1325.    double alpha[2] = {0, 0};
  1326.    double beta[2] = {-1, 0};
  1327.    double A[] = { -0.141, 0.275, 0.717, 0.775, -0.701, -0.689, -0.883, -0.077 };
  1328.    int lda = 4;
  1329.    double B[] = { -0.526, -0.437, 0.133, -0.209, -0.83, 0.328, 0.916, -0.337, 0.762, -0.664, -0.566, 0.955, 0.168, 0.488, -0.172, -0.535 };
  1330.    int ldb = 2;
  1331.    double C[] = { -0.88, 0.945, 0.416, 0.99 };
  1332.    int ldc = 1;
  1333.    double C_expected[] = { 0.88, -0.945, -0.416, -0.99 };
  1334.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1335.    {
  1336.      int i;
  1337.      for (i = 0; i < 2; i++) {
  1338.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1515) real");
  1339.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1515) imag");
  1340.      };
  1341.    };
  1342.   };
  1343.  
  1344.  
  1345.   {
  1346.    int order = 101;
  1347.    int transA = 113;
  1348.    int transB = 113;
  1349.    int M = 1;
  1350.    int N = 2;
  1351.    int K = 4;
  1352.    double alpha[2] = {-0.3, 0.1};
  1353.    double beta[2] = {0, 0.1};
  1354.    double A[] = { -0.534, -0.013, -0.258, -0.31, -0.211, -0.883, -0.89, -0.499 };
  1355.    int lda = 1;
  1356.    double B[] = { -0.185, -0.798, -0.34, 0.716, 0.035, 0.968, -0.26, 0.784, -0.889, -0.344, -0.685, -0.647, -0.764, 0.03, 0.626, -0.989 };
  1357.    int ldb = 4;
  1358.    double C[] = { -0.793, -0.551, 0.182, 0.838 };
  1359.    int ldc = 2;
  1360.    double C_expected[] = { -0.5507177, -0.0286821, 0.2222276, 0.5197398 };
  1361.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1362.    {
  1363.      int i;
  1364.      for (i = 0; i < 2; i++) {
  1365.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1516) real");
  1366.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1516) imag");
  1367.      };
  1368.    };
  1369.   };
  1370.  
  1371.  
  1372.   {
  1373.    int order = 102;
  1374.    int transA = 113;
  1375.    int transB = 113;
  1376.    int M = 1;
  1377.    int N = 2;
  1378.    int K = 4;
  1379.    double alpha[2] = {-0.3, 0.1};
  1380.    double beta[2] = {0, 0.1};
  1381.    double A[] = { 0.575, -0.128, -0.702, 0.758, 0.383, -0.914, 0.157, 0.368 };
  1382.    int lda = 4;
  1383.    double B[] = { 0.572, -0.841, 0.223, -0.334, -0.823, -0.84, 0.671, -0.871, 0.241, 0.927, -0.344, 0.281, -0.034, -0.104, 0.587, -0.329 };
  1384.    int ldb = 2;
  1385.    double C[] = { -0.612, 0.167, 0.647, 0.447 };
  1386.    int ldc = 1;
  1387.    double C_expected[] = { -0.7876717, 0.0341179, -0.0800018, 0.5717566 };
  1388.    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  1389.    {
  1390.      int i;
  1391.      for (i = 0; i < 2; i++) {
  1392.        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1517) real");
  1393.        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1517) imag");
  1394.      };
  1395.    };
  1396.   };
  1397.  
  1398.  
  1399. }
  1400.